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We discuss the quenched limit of lattice QCD at non-zero baryon number den- 
sity. We find evidence for a mixed phase that becomes broader with increasing 
baryon number. Although the action is explicitly Z(3) symmetric the Polyakov 
loop expectation value becomes non-zero already in the low temperature phase. 
It indicates that the heavy quark potential stays finite at large distances, i.e. the 
string between static quarks breaks at non-zero baryon number density already in 
the hadronic phase. This behaviour is validated by calculating the heavy quark 
potential using Polyakov loop correlations. 

1 Introduction 

A quantitative analysis of QCD at non-zero density is important for our un- 
derstanding of the behaviour of dense matter as it is created in heavy ion 
collisions and exists in the cosmological context. While the QCD phase di- 
agram is well known for vanishing baryon density from lattice QCD, for the 
region of non-zero density only qualitative features can be understood in terms 
of models (bags, percolation, strings, ...) and approximations (resonance gas, 
perturbation theory, instanton liquid, ...). Notu-zjero baryon number is usually 
introduced by a non-zero chemical potential /Jim. This leads to a break down 
of the probabilistic interpretation of the path integral representation of the 
QCD partition function as the fermion determinant gets complexH. 
The static limit of QCD at non-zero chemical potential fj, has been formulated 
by Bender et al.ta and Blum et al.B. It seems that in this case the first order 
deconfmement transition of the SU(3) gauge theory turns into a crossover 
for arbitrarily small, non-zero values of the chemical potentials. Rather than 
introducing a non-vanishing chemical potential, i.e. formulate QCD at non- 
vanishing baryon number density in the grand canonical ensemble, one may 
go over to a canonical formulation of the thermodynamics and fix directly 
the baryon naimbcrtl. This is achieved by introducing an imaginary chem- 
ical potentialBQ in the grand canonical partition function and performing a 
Fourier integration to project onto the canonical partition function for a given 
sector of fixed baryon numberl 
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Figure 1. (sgn(Re/g))|| for B = 6 and 12 and lattices of size 8 3 X 2 and 10 3 X 2. 



This leads to a well defined quenched theoryel at fixed baryon number, which 
describes the thermodynamics of gluons in the background of a non-zero num- 
ber of static quark sources, suitably arranged to obey Fermi statistics. 



2 Simulation of quenched QCD with non-zero baryon number 

For any fixed value of the baryon number we can write the quenched partition 
function as 



Z(B,T,V) = J J[d/ •,,,//,.( 



(2) 



where the constraint on the baryon number is encoded in the func- 
tion fg which is a function of Polyakov loops and B counts the number of 
quarks, i.e. B/3 is the baryon number. For B — 3 fg is, for instance, given by 



/ b=3 = (2k) 3 ^(^[L 1)0 ] 3 

+y 2 (8[i 1) o][i 2 ,o]-4[Li,o][io,i]) 
+V(12+^[L 3 ,o]-2[L ll i])) 
with [L itj ] = V- 1 Y, S { r ^L s ) i {TtL%y . 



(3) 
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Figure 2. Polyakov loop expectation value (left) and Polyakov loop susceptibility (right) 
for different values of B and lattices with spatial extend = 8 and 10. 



For a more detailed description of fg see Engels et al.i. Sq is the glu- 
onic action, which is Z(3) symmetric. The partition function Z(B,T,V), is 
non-zero only if B is a multiple of 3, because fg is invariant under Z(3) 
transformations only if B is a multiple of 3. In general it changes by a factor 
e 27riB/3 un( j er a global Z(3) transformation of timelike link variables. 
fs is still a complex function, but upon integration over the gauge fields the 
imaginary part of the partition function vanishes. The remaining sign proh- 
lem can be handled by including the sign in the calculation of observablesu. 
We have performed simulations for the one flavour case (rif — 1) using the 
partition function 

Z {l (B,T,V)= fl[dU x ^\Ref B \e- SG . (4) 

J x,u 

Expectation values of an observable O will be calculated according to 

(O) = (0 - Sgn(R ^ )) ". (5) 
(sgn(Re/ B ))|| 

The simulations are performed on N 3 xN T lattices with N a — 8, 10 and N T = 2 
using the standard Wilson action with quark number values of B = 6 and 12. 
When increasing the gauge coupling (3 the temperature is increased while 
ub/T 3 = ^B(j^-) 3 is kept fixed. Close to T c a simulation on a 8 3 x 2 lattice 
with B — 12 corresponds to tib — 0.15//m 3 , i.e. approximately nuclear 
matter density. Figure 1 shows the average sign (sgn(/s))|| as a function of 
the coupling /3. For large values of the temperature the sign is almost always 
positive, but also for the smallest temperature in our analysis the sign can be 
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Figure 3. Heavy quark potential for T ~ 0.86T C and B = and 6. 



well determined. It depends on the spatial volume N% but varies little with 
B. 

In the quenched limit at vanishing baryon number the Polyakov loop ex- 
pectation (L) value is an order parameter for the deconfinement transition in 
the infinite volume limit. For a first order phase transition, (L) changes dis- 
continuously at T c . In Figure 2 we see a clear signal for a first order transition 
for the B — case, while for all B > the transition is continuous. The tran- 
sition region is shifted towards smaller values and it broadens with increasing 
B. This is the expected behaviour for a canonical calculation. By changing 
the gauge coupling (3 we vary the lattice cut-off and through this also the 
baryon number density continuously. At fixed baryon number we therefore 
follow a simulation path that traverses the mixed phase continuously. 
If there is a mixed phase in the sense that there exist singularities in thermo- 
dynamic observables when entering and leaving this mixed phase, it could be 
reflected in a discontinuous change of the slope of the Polyakov loop expec- 
tation value. We analyzed this by calculating the conventional Polyakov loop 
susceptibility \l (Figure 2). This response function reflects the existence of 
a transition region that becomes broader with increasing ub, but does not 
show indications for a discontinuity. 

The Polyakov loop expectation value becomes non-zero already in the 
low temperature phase. This indicates that the heavy quark potential stays 
finite at large distances. We validate this by calculating the potential using 
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Figure 4. Plaquette expectation value (left) and the negative of the free energy density 
(right). On the right plot we only show tick marks at critical couplings on lattices with 
extend N T = 2, 3, 4, 6 and 8, which corresponds here to temperatures T/T c = 1, 1.5, 2, 3, 4 



Polyakov loop correlations. The potentials for T ~ 0.86T C at B = and 6 on 
a 16 3 x 4 lattice are plotted in Figure 3. For zero baryon number it shows 
the usual behaviour for the quenched case. The potential is linearly rising 
for large distances. For B — 6 the potential stays finite at large distances. 
The static quark anti-quark sources used to probe the heavy quark potential 
can recombine with the already present static quarks. This leads-,to string 



The plaquette expectation value (Figure 4) shows a similar behaviour 
as the Polyakov loop. With increasing ns the transition region broadens. 
Through an integration over differences of plaquette expectation values for 
finite temperature and zero temperature, the free energy density can be cal- 
culated. As the area between these data increases with increasing jib, the 
free energy density decreases at fixed temperature with increasing B. 



3 Conclusions 

We have analyzed the quenched limit of QCD at non-zero baryon number. 
Although a sign problem remains in this theory, it can be handled quite well 
numerically. We find indications for a mixed phase, which broadens with 
increasing baryon number density and is shifted towards smaller tempera- 
tures. We also see evidence that the heavy quark potential stays finite for 
large distances already in the hadronic phase. String breaking starts at short 
distances. 



baryon: submitted to World Scientific on February 1, 2008 5 




Acknowledgments 

This work was partly supported by the TMR network Finite Temperature 

Phase Transitions in Particle Physics, EU contract no. ERBFMRX-CT97- 
0122 and the DFG through grant no. KA 1198/4-1. 



References 



1. P. Hasenfratz and F. Karsch, Phys. Lett. 125B (1983) 308. 

2. J. Kogut, M. Matsuoka, M. Stone, H.W. Wyld, J.H. Shenker, J. 
Shigemitsu and D.K. Sinclair, Nucl. Phys. B225 [FS9] (1983) 93. 

3. For a recent review see : I.M. Barbour, S.E. Morrison, E.G. Klepfish, J.B. 
Kogut and M.P. Lombardo, Nucl. Phys. B (Proc. Suppl.) 60 (1998) 220. 

4. I. Bender, T. Hashimoto, F. Karsch, V. Linke, A. Nakamura, M. Plewnia, 
I.O. Stamatescu and W. Wetzel, Nucl. Phys. B (Proc. Suppl.) 26 (1992) 
323. 

5. T.C. Blum, J.E. Hetrick and D. Toussaint, Phys. Rev. Lett. 76 (1996) 
1019. 

6. D.E. Miller and K. Redlich, Phys. Rev. D35 (1987) 2524. 

7. A. Roberge and N. Weiss, Nucl. Phys. B275 [FS17] (19 86) 734. 

8. J. Engels, O. Kaczmarek, F. Karsch and E. Laermann, iicp-lat /9903030[ 

9. B. Berg, J. Engels, E. Kchl, B. Waltl and H. Satz, Z. Phys. C31 (1986) 
167. 

10. C. DeTar, O. Kaczmarek, F. Karsch and E. Laermann, Phys. Rev. D59 
(1999) 031501. 



baryon: submitted to World Scientific on February 1, 2008 6 



